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Abstract 

The decays and are studied with the ATLAS detector at the 

LHC using a dataset corresponding to integrated luminosities of 4.9 fb”* and 20.6 fb”* 
of pp collisions collected at centre-of-mass energies ^[s - 1 TeV and 8 TeV, respect¬ 
ively. Signal candidates are identified through Jjip p^p~ and (pn'^{yln^) de¬ 

cays. With a two-dimensional likelihood fit involving the B^ reconstructed invariant mass 
and an angle between the and candidate momenta in the muon pair rest frame, the 
yields of B^ f^e transverse polarisation fraction in 

B^ decay are measured. The transverse polarisation fraction is determined to 

be r++(B^ ^ = 0.38 + 0.23 + 0.07, and the derived ratio of the 

branching fractions of the two modes is 3b*^ji,i,dp = 2.8;|;Qg + 0.3, where the 
first error is statistical and the second is systematic. Finally, a sample of B^ decays 

is used to derive the ratios of branching fractions 3b*^ji,i/d* lfBB+->ji^n+ = 3.8 +1.1 + 0.4+0.2 
and Sb+->ji^d’PISb*->jii/, 7 t+ = 10.4 + 3.1 + 1.5 + 0.6, where the third error corresponds to the 
uncertainty of the branching fraction of <p{K'^K~)n'^ decay. The available theoretical 

predictions are generally consistent with the measurement. 
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1 Introduction 


The meson^ is the only known weakly decaying particle consisting of two heavy quarks. The ground 
be state was first observed by CDF [1] via its semileptonic decay ^ Jlxpi^vt. An excited be state 
has been observed recently by ATLAS [2] using the B* decay mode . The presence of two 

heavy quarks, each of which can decay weakly, affects theoretical calculations of the decay properties of 
the meson. In the case of ^ ^ ccs processes, decays to charmonium and a or a meson are 
predicted to occur via colour-suppressed and colour-favoured spectator diagrams as well as via the weak 
annihilation diagram (see Fig. 1). The latter, in contrast to decays of other B mesons, is not Cabibbo- 
suppressed and can contribute significantly to the decay amplitudes. The decay properties are addressed 
in various theoretical calculations [3-9] and can also be compared to the analogous properties in the 
lighter B meson systems such as B|J ^ or B^ The decays Bj ^ Jjij/D* and 

B^ ^ Jlif/D*'*', which have been observed recently by the LHCb experiment [10], provide a means to test 
these theoretical predictions. 



Figure 1: Feynman diagrams for B^ —> JI\j/D^P^ decays: (a) colour-favoured spectator, (b) colour-suppressed 
spectator, and (c) annihilation topology. 

This paper presents a measurement of the branching fractions of B^ ^ and B^ ^ JjipD** decays, 

normalised to that of B^ ^ decay, and polarisation in B^ ^ Jfif/D*'^ decay performed with the 

ATLAS detector [11]. The meson is reconstructed via the (pn''' decay with the ^ meson decaying 

into a pair of charged kaons. The D*^ meson decays into a DJ meson and a soft photon or Detecting 
such soft neutral particles is very challenging, thus no attempt to reconstruct them is made in the analysis. 
The 7/i/^ meson is reconstructed via its decay into a muon pair. 

The measurement presented in this paper allows an independent verification of the results of Ref. [10] 
with comparable statistical and systematic uncertainties. The following ratios are measured: = 

and!/?£)•+//>+ = I^Bp^JiiijD+^ where 

!Bb+->x denotes the branching fraction of the B^ ^ X decay. The decay B^ —> Jft//D*^ is a transition of 
a pseudoscalar meson into a pair of vector states and is thus described by the three helicity amplitudes, 
A++, A _, and Aqo, where the subscripts correspond to the helicities of //i/^ and D*^ mesons. The contri¬ 
bution of the A++ and A_amplitudes, referred to as the A++ component, corresponds to the Jlip and D*'*' 

transverse polarisation. The fraction of transverse polarisation, F++/F = F++(B^ ^ 7/(AD*''')/F(B^ ^ 


* Charge conjugate states are implied throughout the paper unless otherwise stated. 
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Jlij/D**), is also measured. From a naive prediction by spin counting, one would expect this fraction to 
be 2/3, while calculations [8, 9] predict values of 0.41-0.48. 

This analysis is based on a combined sample of pp collision data collected by the ATLAS experiment at 
the LHC at centre-of-mass energies ^[s - 1 TeV and 8 TeV corresponding to integrated luminosities of 
4.9 fb“^ and 20.6 fb“^ respectively. 

2 The ATLAS detector, trigger selection and Monte Carlo samples 

ATLAS is a general-purpose detector consisting of several subsystems including the inner detector (ID), 
calorimeters and the muon spectrometer (MS). Muon reconstruction makes use of both the ID and the MS. 
The ID comprises three types of detectors: a silicon pixel detector, a silicon microstrip semiconductor 
tracker (SCT) and a transition radiation tracker. The ID provides a pseudorapidity^ coverage up to \t]\ = 
2.5. Muons pass through the calorimeters and reach the MS if their transverse momentum, px, is above 
approximately 3GeV.^ Muon candidates are formed either from a stand-alone MS track matched to an 
ID track or, in case the MS stand-alone track is not reconstructed, from an ID track extrapolated to the 
MS and matched to track segments in the MS. Candidates of the latter type are referred to as segment- 
tagged muons while the former are called combined muons. Muon track parameters are taken from the 
ID measurement alone in this analysis, since the precision of the measured track parameters for muons in 
the pt range of interest is dominated by the ID track reconstruction. 

The ATLAS trigger system consists of a hardware-based Level-1 trigger and a two-stage High Level 
Trigger (HLT). At Level-1, the muon trigger uses dedicated MS chambers to search for patterns of hits 
satisfying different px thresholds. The region-of-interest around these hit patterns then serves as a seed for 
the HLT muon reconstruction, in which dedicated algorithms are used to incorporate information from 
both the MS and the ID, achieving a position and momentum resolution close to that provided by the 
offline muon reconsfrucfion. Muons are efficienfly friggered in fhe pseudorapidify range \r]\ < 2.4. 

Triggers based on single-muon, dimuon, and fhree-muon signafures are used fo selecf Jlip ^ p*p~ decays 
for fhe analysis. The fhird muon can be produced in fhe signal evenfs in semilepfonic decays of fhe 
fwo ofher heavy-flavour hadrons. The majorify of evenfs are collecfed by dimuon friggers requiring a 
verfex of fwo opposifely charged muons wifh an invarianf mass befween 2.5 GeV and 4.3 GeV. During fhe 
dafa faking, fhe px threshold for muons in these triggers was either 4 GeV or 6 GeV. Single-muon triggers 
additionally increase the acceptance for asymmetric 7/i/^ decays where one muon has px < 4 GeV. Finally, 
three-muon triggers had a px threshold of 4 GeV, thus enhancing the acceptance during the periods of high 
luminosity when the px threshold for at least one muon in the dimuon triggers was 6 GeV. 

Monte Carlo (MC) simulation is used for the event selection criteria optimisation and the calculation 
of the acceptance for the considered decay modes. The MC samples of the B^ decays were gener¬ 
ated with Pythia 6.4 [12] along with a dedicated extension for the B^ production based on calculations 
from Refs. [13-16]. The decays of B^ are then simulated with EvtGen [17]. The generated events were 
passed through a full simulation of the detector using the ATLAS simulation framework [18] based on 
Geant 4 [19, 20] and processed with the same reconstruction algorithms as were used for the data. 

^ ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the i/-axis points 
upward. Cylindrical coordinates (r, (p) are used in the transverse plane, 0 being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 0 as ?/ = - In tan(6/2). 

^ Using a system of units with c = 1 is implied throughout the paper. 
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3 Reconstruction and event selection 


The Jlip candidates are reconstructed from pairs of oppositely charged muons. At least one of the two 
muons is required to be a combined muon. Each pair is fitted to a common vertex [21]. The quality of 
the vertex fit must satisfy ;t'^/ndf <15, where the ndf stands for the number of degrees of freedom. The 
candidates in the invariant mass window 2800 MeV < < 3400MeV are retained. 

For the D* cl){K^K~)n* reconstruction, tracks of particles with opposite charges are assigned kaon 
mass hypotheses and combined in pairs to form (p candidates. An additional track is assigned a pion 
mass and combined with the (p candidate to form a Dj candidate. To ensure good momentum resolution, 
all three tracks are required to have at least two hits in the silicon pixel detector and at least six hits 
in the SCT. Only three-track combinations successfully fitted to a common vertex with /ndf < 8 are 
kept. The (/> candidate invariant mass, m{K^K~), and the D+ candidate invariant mass, m{K^K~n^), are 
calculated using the track momenta refitted to the common vertex. Only candidates with m{K'^K~) within 
±7 MeV around the (p mass, = 1019.461 MeV [22], and with 1930 MeV < m{K*K~n*) < 2010 MeV 
are retained. 

The candidates are built by combining the five tracks of the Jjip and DJ candidates. The 

Jlip meson decays instantly at the same point as the B* does (secondary vertex) while the Dj lives long 
enough to form a displaced tertiary vertex. Therefore the five-track combinations are refitted assuming 
this cascade topology [21]. The invariant mass of the muon pair is constrained to the Jjtp mass, mj/^ = 
3096.916 MeV [22]. The three daughter tracks are constrained to a tertiary vertex and their invariant 
mass is fixed to the mass of D^, mD+ = 1968.30MeV [22]. The combined momentum of the refitted 
decay tracks is constrained to point to the dimuon vertex. The quality of the cascade fit must satisfy 
< 3. 

The B^ meson is reconstructed within the kinematic range p^iB^) > 15 GeV and 177 ( 5^)1 < 2.0, where the 
detector acceptance is high and depends weakly on pt(®c) )• 

The refitted tracks of the daughter hadrons are required to have I 77 I < 2.5 and px > 1 GeV, while the 
muons must have |? 7 | < 2.3 and px > 3 GeV. To further discriminate the sample of D* candidates from a 
large combinatorial background, the following requirements are applied: 

• cos0*(7r) < 0.8, where 6*{n) is the angle between the pion momentum in the rest frame 

and the K"^K~n^ combined momentum in the laboratory frame; 

• I cos^ d'{K)\ > 0.15, where 6'{K) is the angle between one of the kaons and the pion in the K*K~ 
rest frame. The decay of the pseudoscalar DJ meson to the cp (vector) plus n (pseudoscalar) final 
state results in an alignment of the spin of the ^ meson perpendicularly to the direction of motion 
of the (p relative to DJ. Consequently, the distribution of cos6'{K) follows a cos^6'{K) shape, 
implying a uniform distribution for cos^ 0'{K). In contrast, the cos6'{K) distribution of the com¬ 
binatorial background is uniform and its cos^ 0'{K) distribution peaks at zero. The cut suppresses 
the background significantly while reducing the signal by 15%. 

The B^ candidate is required to point back to a primary vertex such that (/^^(B/') < 0.1 mm and 
Zq'^(Bc) sin0(B]() < 0.5 mm, where and are respectively the transverse and longitudinal im¬ 
pact parameters with respect to the primary vertex. All primary vertices in the event are considered. If 
there is more than one primary vertex satisfying these requirements (~ 0.5% events both in data and MC 
simulation), the one with the largest sum of squared transverse momenta of the tracks originating from it 
is chosen. 
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The transverse decay length"^ of the candidate is required to satisfy Lxy{B"^) >0.1 mm. The transverse 
decay length of the D+ measured from the B"^ vertex must be L^yiD'D > 0.15 mm. In order to remove 
fake candidates, both Lxy{B'^) and Lxy{D^) are required not to exceed 10 mm. 

Taking into account the characteristic hard fragmentation of Zj-quarks, a requirement pt{B*)I Yj PT(trk) > 
0.1 is applied, where the sum in the denominator is taken over all tracks originating from the primary 
vertex (tracks of the B^ candidate are included in the sum if they are associated with the primary vertex). 
The requirement reduces a sizeable fraction of combinatorial background while having almost no effect 
on the signal. 

The following angular selection requirements are introduced to further suppress the combinatorial back¬ 
ground: 

• cos0*(Dj) > -0.8, where 0*(DJ) is the angle between the DJ candidate momentum in the rest 

frame of the candidate, and the B^ candidate line of flight in the laboratory frame. The dis¬ 
tribution of cos0*(Dj) is uniform for the decays of pseudoscalar meson before any kinematic 
selection while it tends to increase for negative values of cos for the background. 

• cos6'{n) > -0.8, where 6'{n) is the angle between the Jj^ candidate momentum and the pion 
momentum in the K^K~n^ rest frame. Its distribution is nearly uniform for the signal processes 
but peaks towards -1 for the background. 

Distributions of these two variables after applying all other selection requirements described in this sec¬ 
tion are shown in Fig. 2. They are shown for the simulated signal samples, as well as for sidebands of 
the mass spectrum in data, defined as the regions 5640 MeV < m{Jlij/D*) < 5900 MeV (left sideband) 
and 6360 MeV < m{JliJ/D^) < 6760 MeV (right sideband). A dip in the cos0'(7r) distribution for the 
B^ signal is caused by rejection of B® ^ candidates discussed below. 




Figure 2: Distributions of (a) cos and (b) cos 9'{n), where 9*(D^) and 9'(7t) are two angular variables defined 

in Sect. 3. The distributions are shown for data sidebands (black dots) and MC simulation of —> Jl^fiD^ signal 
(red solid line) and Aqo (green dotted line) and A±± (blue dashed line) components of JfiJ/DY signal. The 

distributions are obtained after applying all selection criteria except the ones on the plotted variable. The MC 
distributions are normalised to data. 


The transverse decay length of a particle is defined as the transverse distance between the production (primary) vertex and the 
particle decay (secondary) vertex projected along its transverse momentum. 
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Various possible contributions of partially reconstructed B JjijjX decays were studied. The only signi¬ 

ficant one was found from the ^ decay process. This contribution arises when the combination 
of the tracks from a true ^ J^~) decay with a fifth random track results in a fake 
—> JltJ/(jj.*ii~)D^(K^K~7T'*') candidate. For each reconstructed B^ candidate, an additional vertex fit 
is performed. The two muon tracks and the two kaon tracks are fitted to a common vertex, where the 
kaon tracks are assumed to be from (p K^K~ and the muon pair is constrained to have the nominal Jjifj 
mass. The mass of the candidate, m{ii^K~), is then calculated from the refitted track parameters. 
Candidates with 5340 MeV < jTK~) < 5400MeV are rejected. This requirement suppresses the 
bulk of the events while rejecting only ~ 4% of the signal. 

After applying the selection requirements described above, 1547 JlipD^ candidates are selected in the 
mass range 5640-6760 MeV. 

4 J candidate fit 

The mass dishibution of the selected B^ candidates is shown in Fig. 3. The peak near 

the B^ mass, ms+ = 6275.6MeV [22], is athibuted to the signal of B^ decay while a wider 

shucture between 5900 MeV and 6200 MeV corresponds to B^ JlipD** with subsequent D** D'^y 
or D*^ decays where the neutral particle is not reconstructed. 



m(J/i|/Dp [MeV] 

Figure 3: The mass distribution for the selected candidates. The red solid line represents the projection of 

the fit to the model described in the text. The contribution of the decay is shown with the magenta 

long-dashed line; the brown dash-dot and green dotted lines show the Aqo and component 

contributions, respectively; the blue dashed line shows the background model. The uncertainties of the listed fit 
result values are statistical only. 

Mass distributions of the JItp and candidates corresponding to the JliJ/D^ mass region of the ob¬ 
served B^ signals are shown in Fig. 4. To obtain these plots, the B^ candidates are built 

without the mass constraints in the cascade fit, with the mass of the candidate calculated as nr(7/i/rDj) = 
jTyr)+mji^-m{K^K~n^)+mD+, where andm£)+ are the nominal masses of the 
respective particles. The mass of the Bp candidate is required to be 5900 MeV < m{JlipD^) < 6400 MeV 
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while the mass windows for the corresponding intermediate resonances are widened to the plotting ranges. 
The Jjip and DJ mass distributions are fitted with a sum of an exponential function describing the back¬ 
ground and a modified Gaussian function [23,24] describing the corresponding signal peak. The modified 
Gaussian function is defined as 


Gauss“°‘^ ~ exp 


1 + 


l+.t/2 


( 1 ) 


where x = |mo - m|/cr with the mean mass mo and width cr being free parameters. The fitted masses 
of JI4/ (3095.1 ± 2.4MeV) and (1969.0 ± 4.1 MeV) agree with their nominal masses, the widths are 
consistent with those in the simulated samples, and the signal yields are found to h&Njj^ = 568 ± 28 and 
Nd*- - 175 ± 36. 



Figure 4; Mass distribution of the (a) Jj^ and (b) candidates after the full selection (without 

mass constraints in the cascade fit) in the mass window of the candidate 5900 MeV < m(Jlil/D^) < 6400 MeV. 
The spectra are fitted with a sum of an exponential and a modified Gaussian function. The uncertainties of the 
shown 7/i/f and yields are statistical only. 


The information about the helicity in JI decay is encoded both in the mass distribution of 

the JI system and in the distribution of the helicity angle, which is defined in the rest frame 

of the muon pair as the angle between the /r'*' and the Dj candidate momenta. Thus, a two-dimensional 
extended unbinned maximum-likelihood fit of the m{Jft//D^) and | cos O'Qi'*')] distributions is performed. 

The A++ and A_helicity amplitude contributions are described by the same mass and angular shapes 

because of the parity symmetry of the 7/(/^ and decays. This is confirmed by the MC simulation. 
Thus these components are treated together as the A++ component, while the shape of the Aqo component 
is different and is therefore treated separately. A simultaneous fit to the mass and angular distributions 
significantly improves the sensitivity to the contributions of the helicity amplitudes in —> Jlt^D*'^ 

decay with respect to a one-dimensional mass fit. 

Four two-dimensional probability density functions (PDFs) are defined to describe the B^ 
signal, the A++ and Aqo components of the B^ —> JftJ/D*^ signal, and the background. The signal PDFs 
are factorised into mass and angular components. The effect of correlations between their mass and 
angular shapes is found to be small and is accounted for as a systematic uncertainty. 

The mass distribution of the B^ signal is described by a modified Gaussian function. For the 

B^ JlipD*'^ signal components, the mass shape templates obtained from the simulation with the kernel 
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estimation technique [25] are used. The branching fractions of D** and D*'*' D'^y decays for 

the simulation are set to the world average values [22]. The position of the templates along the mass axis is 
varied in the fit simultaneously with the position of the JItpD'l signal peak. The background mass 

shape is described with a two-parameter exponential function, exp [a • + b ■ 

To describe the | cos 0'(ju''‘)| shapes, templates from the kernel estimation are used. The templates for 
the signal angular PDFs are extracted from the simulated samples. Although their shapes are calculable 
analytically, using the templates allows the fit to account for detector effects. The background angular 
description is based on the | cos 0'(jj.*)\ shape of the candidates in the sidebands of JliJ/D^ mass spectra. 
Two templates are produced from the angular distributions of the candidates in the left and right mass 
sidebands as defined in Sect. 3. The angular PDF for the background is defined as a conditional PDF of 
I cos 6'Qi''')\ given the per-candidate m{JlipD^). For the candidates in the lower half of the left sideband 
(5640-5770 MeV), the template from the left sideband is used. Similarly, the template from the right 
sideband is used for the upper half of the right sideband (6560-6760 MeV). For the candidates in the 
middle part of the mass spectrum (5770-6560 MeV), a linear interpolation between the two templates is 
used. 

The fit has seven free parameters: the mass of the meson, the relative contribution of the 

A++ component to the total B^ Jjif/D** decay rate in the selected sample, /+±; the two parameters 
of the exponential background; the yields of the two signal modes, Ng+^jjij^Dy, and the 

background yield. The width of the modified Gaussian function, o-b+^j/i/,d+-: is fixed fo the value obtained 
from the fit to the simulated signal, crs+^j/if,D+ = 9.95 MeV. Leaving this parameter free in the data fit 
results in the value 7.9 ± 3.0 MeV, consistent with the simulation in the range of statistical uncertainty. 

It was checked that the fit procedure provides unbiased values and correct statistical uncertainties for the 
extracted parameters using pseudo-experiments. The values of the relevant parameters obtained from 
the fit are given in Table 1. The fitted B^ mass agrees with the world average value [22]. The mass 
and angular projections of the fit on the selected candidate dataset are also shown in Figs. 3 

and 5(a), respectively. In order to illustrate the effect of the angular part of the fit in separating the 
helicity amplitudes, the | cos O'iju''')] projection for the subset of candidates with the masses 5950 MeV < 
< 6250 MeV corresponding to the region of the observed B^ JjiJ/D** signal is shown in 

Fig. 5(b). 


Table 1: Parameters of the JI signals obtained with the unbinned extended maximum-likelihood fit. 

The width parameter of the modified Gaussian function is fixed to the MC value. Only statistical uncertainties are 
shown. No acceptance corrections are applied to the signal yields. 


Parameter 

Value 

1 ' 1 /D* [MeV] 

6279.9 ± 3.5 


36 ± 10 

NBt^JI<pDT 

95 ±27 

f±± 

0.37 ± 0.22 


The statistical significance for the observed B^ signal estimated from toy MC studies is 4.9 standard 
deviations. 









Figure 5; The projection of the likelihood fit on the variable | cos where the helicity angle is the angle 

between the and candidate momenta in the rest frame of the muon pair from 7/i/f decay, for (a) the full 
selected candidate dataset and (b) a subset of the candidates in a mass range 5950 MeV < m(Jftl/D^) < 

6250 MeV corresponding to the observed signal of decay. The red solid line represents the full ht 

projection. The contribution of the decay is shown with the magenta long-dashed line (it is not drawn 

in (b) because this contribution vanishes in that mass range); the brown dash-dot and green dotted lines show the 
B^ Aoo and A±± component contributions, respectively; the blue dashed line shows the background 

model. 

S J 1^71'^ candidate reconstruction and fit 

JjijjT:* candidates are reconstructed by fitting a common vertex of a pion candidate track and the 
two muons from a Jjij/ candidate, selected as described in Sect. 3. For the pion candidate, tracks identified 
as muons are vetoed in order to suppress the substantial background from —> Jlipii^VnX decays. The 
invariant mass of the two muons in the vertex fit is constrained to the Jji}/ nominal mass. The quality 
of the fit must satisfy ;y^/ndf < 3. The following selection requirements applied to the JjijjT:* 

candidates are analogous to those for B^ JIipD* candidates described in Sect. 3: the candidates must 
be within the kinematic range > 15GeV, \r]iB'^)\ < 2.0; the refitted values of the transverse 

momenta and pseudorapidities of the muons are required to satisfy > 3GeV, \T]ip.-)\ < 2.3; the 

same requirements on pointing to the primary vertex and the ratio )/FT(trk) are applied. The 

refitted pion track kinematics must satisfy > 5 GeV and |?7(7r''‘)| < 2.5. The transverse decay length 

is required to be Lxy{B'^) > 0.2 mm, and not to exceed 10 mm. 

To further suppress combinatorial background, the following selection is applied: 

• cos > -0.8, where 6*{n) is the angle between the pion momentum in the rest frame 

and the B^ candidate line of flight in laboratory frame. This angular variable behaviour for the 
signal and the background is the same as that of cos 0*(Dj) used for Jlij/D* candidates selection. 

• I cos 9'{p*)\ < 0.8, where is the angle between the and momenta in the muon pair rest 

frame. The signal distribution follows a sin^ d'ip*) shape, while the background is flat. 

After applying the above-mentioned requirements, 38542 candidates are selected in the mass range 

5640-6760 MeV. Figure 6 shows the mass distribution of the selected candidates. An extended unbinned 
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maximum-likelihood fit of the mass speetrum is performed to evaluate the signal yield. The 

signal contribution is described with the modihed Gaussian function while an exponential function is used 
for the background. The mass, the width of the modihed Gaussian function, crB+^ji,ijn+^ 

the yields of the signal, Nb+^ji^j^+, and the background, and the slope of the exponential background are 
free parameters of the ht. The ht results are summarised in Table 2, and the ht projection is also shown in 
Fig. 6. The extracted B^ mass value is consistent with the world average [22], and the signal peak width 
agrees with the simulation (37.4 MeV). 


Table 2: Signal parameters of the mass distribution obtained with the unbinned extended maximum- 

likelihood fit. Only statistical uncertainties are shown. No acceptance corrections are applied to the signal yields. 


Parameter 

Value 

[MeV] 

6279.9 ± 3.9 

o'b+^jI4/7t+ [MeV] 

33.9 ± 4.2 

NB*^JI^n+ 

1140 ± 120 



Figure 6; The mass disttibution for the selected candidates. The red solid line represents the result of 

the ht to the model described in the text. The brown dotted and blue dashed lines show the signal and background 
component projections, respectively. The uncertainty of the shown signal yield is statistical only. 


6 Branching fractions and polarisation measurement 

The ratios of the branching fractions 7?/)+/;r+ and are calculated as 

/<^(»)+, + — —— - — —- X - X -, (Zj 

where JiB+-,x and Nb+-,x are the total acceptance and the yield of the corresponding mode. For Sj)+^(f,{K+K-) 7 T+^ 
the CLEG measurement [26] of the partial Dj ^ K^K~7t^ branching fractions, with a kaon-pair mass 
within various intervals around the nominal (p meson mass, is used. An interpolation between the partial 
branching fractions, measured for ±5 MeV and ± 10 MeV intervals, using a relativistic Breit-Wigner shape 
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of the resonance yields the value (1.85 + 0.11)% for the ±7MeV interval which is used in the analysis. 
The effect of admixture of other DJ decay modes with {K"^K~n^) final state which are not present in the 
MC simulation is studied separately and accounted for as a systematic uncertainty. 

The acceptance for the Jjif/D** decay mode is different for the A±± and Aqo components, thus the 
full acceptance for the mode is 

^B+^JI4,DA = - — - + - — -) > (3) 

\^Bt^Jlil>DA,A±± ^Bt^Jl^lyDAAool 

where the subscripts indicate the helicity state and /++ is the value extracted from the fit (Table 1). The 
acceptances are determined from the simulation and shown in Table 3. 


Table 3; The acceptance for all decay modes studied. Only uncertainties due to MC statistics are shown. 


Mode 

[%] 

B+ ^ 7/tA7r+ 

4.106 ±0.056 

^ 7 / 1 AD+ 

1.849 ±0.034 


1.829 ±0.053 

B+ ^ A±± 

1.712 ±0.035 


The ratio ‘Rd*+id^, is calculated as 


'f<DAIDt 


^B+^JlijjDA 


Nb+^JIiPD*+ ^Bt^JlijjD* 

X j 

Nb+^JI^D+ ^Bt^JlijjDA 


(4) 


where the ratio of the yields its uncertainty is extracted from the fit as a 

parameter in order to account for correlations between the yields. 


The fraction of the A++ component contribution in B* —> JlipD*'^ decay is calculated from the f±± value 
quoted in Table 1 by applying a correction to account for the different acceptances for the two component 
contributions: 


r±±/r = /±±x 


^Bt^JI<ljDT 

^B+^JIiliD*+,^±± 


(5) 


7 Systematic uncertainties 


The systematic uncertainties of the measured values are determined by varying the analysis procedure 
and repeating all calculations. Although some sources can have rather large effects on the individual 
decay rate measurements, they largely cancel in the ratios of the branching fractions due to correlation 
between the effects on the different decay modes. The following groups of systematic uncertainties are 
considered. 

The first group of sources of systematic uncertainty relates to possible differences between the data and 
simulation affecting the acceptances for the decay modes. Thus, an effect of the B* production model is 
evaluated by varying the simulated pj and \t]\ spectra while preserving agreement with the data distribu¬ 
tions obtained using the abundant B^ channel. These variations have very similar effects on 
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the acceptances for the different decay modes, thus giving rather moderate estimates of the uncertainties, 
not exceeding 3% in total, on the ratios of branching fractions. The effect of presence of other decay 
modes with {K^K~n^) final state on the calculated acceptances is studied with a separate MC simulation. 
Its conservative estimate yields 0.4% which is assigned as *??£)+ and ’/?£)•+/;j+ uncertainties. An uncer¬ 
tainty on the tracking efficiency is dominated by the uncertainty of the detector material description in 
the MC simulation. Samples generated with distorted geometries and with increased material are used 
to estimate the effect on track reconstruction efficiencies. When propagated to the ratios of branching 
fractions, these estimates give 0.5% uncertainty for /?i-+ tracks in 

modes. Limited knowledge of the and DJ lifetimes leads to an additional systematic 
uncertainty. The simulated proper decay times are varied within one standard deviation from the world 
average values [22] resulting in uncertainties of ~ 1% assigned to ^D*+in+ due to the B^ 

lifetime, and 0.3% due to the lifetime. Removing the requirement on )/ ^ /^T(trk) is found to 
produce no noticeable effect on the measured values. 

The next group of uncertainties originates from the signal extraction procedure. These uncertainties are 
evaluated separately for 7/i/^Dj and JjijjT:* candidate fits. For the former, the following variations of the 
fit model are applied and the difference is treated as a systematic uncertainty: 

• different background mass shape parametrisations (three-parameter exponential, second- and third- 
order polynomials), different fitted mass range (reduced by up to 40MeV from each side independ¬ 
ently); 

• a double Gaussian or double-sided Crystal Ball function [27-29] for > Jlij/D^ signal descrip¬ 
tion; variation of the modified Gaussian width within 10% of the MC simulation value; 

• variation of the smoothness of the B^ JliJ/D** signal mass templates, which is controlled by a 
parameter of the kernel estimation procedure [25] ; 

• similar variation of the smoothness of the B^ —> signal angular templates; 

• variation of the smoothness of the sideband templates used for the background angular PDF con¬ 
struction; different ranges of the sidebands; different sideband interpolation procedure; 

• modelling of the correlation between the mass and angular parts of the signal PDFs. This correla¬ 

tion takes place only at the detector level and manifests itself in degradation of the mass resolution 
for higher values of | cos A dedicated fit model accounting for this effect is used for the data 

fit. The impact on the result is found to be negligible compared to the total uncertainty. 

The first two items give the dominant contributions to the uncertainties of the ratios of branching fractions 
while the transverse polarisation fraction measurement is mostly affected by the background angular 
modelling variations. For the normalisation channel fit model, the similar variations of the background 
and signal mass shape parametrisation are applied. The deviations produced by the variations of the fits 
reach values as high as 10-15% thus making them the dominant sources of systematic uncertainty. 

The branching fractions of D*"*" [22] are varied in simulation within their uncertainties to estimate their 
effect on the measured quantities. Very small uncertainties are obtained for the 'Rd’+/ 7 t+ and 'Rd-*id+^ 
while for F++/F, the estimate is ~ 1%. 

The statistical uncertainties on the acceptance values due to the MC sample sizes are also treated as a 
separate source of systematic uncertainty and estimated to be 2-3%. 
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In order to check for a possible bias from using three-muon triggers, vetoing the meson daughter 
tracks identified as muons is tested and found not to affect the measurement. 

Finally, since Sd+_, 0 (ic+k-) 7 t* enters Eq. (2), its uncertainty, evaluated from Ref. [26] as 5.9%, is propag¬ 
ated to the final values of the relative branching fractions. 

The systematic uncertainties on the measured quantities are summarised in Table 4. 

Table 4; Relative systematic uncertainties on the measured ratios of branching fractions Rd*/ 7 r+, RD‘;*/jr+, Rd'+/d* and 
on the transverse polarisation fraction r±±/r. 


Source 


Uncertainty [%] 

Ru-pln* 

r±±/r 

Simulated pjiB^) spectrum 

0.4 

0.9 

0.5 

0.4 

Simulated spectrum 

1.9 

2.4 

0.6 

0.2 

Other Dj decay modes contribution 

0.4 

0.4 

- 

- 

Tracking efficiency 

0.5 

0.5 

<0.1 

<0.1 

B^ lifetime 

1.2 

1.3 

<0.1 

<0.1 

Dj lifetime 

0.3 

0.3 

<0.1 

<0.1 

B^ signal extraction 

4.4 

10.5 

10.7 

17.4 

B^ signal extraction 

8.5 

8.5 

- 

- 

D*'^ branching fractions 

<0.1 

< 0.1 

<0.1 

1.1 

MC sample sizes 

2.3 

2.4 

2.7 

2.2 

Total 

10.1 

14.0 

11.0 

17.6 

^D+-kI>{K+K-)7t+ 

5.9 

5.9 

- 

- 


8 Results 


The following ratios of the branching fractions are measured: 


^ ^ = 3.8 ± 1.1 (stat.) ± 0.4 (syst.) ± 0.2 (BE), 

= 10-4 ±3.1 (stat.) ± 1.5 (syst.) ± 0.6 (BE), 
= l.Slli (Stat.) ± 0.3 (syst.). 


( 6 ) 

(V) 

( 8 ) 


where the BE uncertainty corresponds to the knowledge of SD+^(f,{K+K-)n+- The relative contribution of 
the A++ component in > JjiJ/D*'*' decay is measured to be 


r±±/r - 0.38 ± 0.23 (stat.) ± 0.07 (syst.) 


(9) 
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These results are compared with those of the LHCb measurement [10] and to the expectations from 
various theoretical calculations in Table 5 and Fig. 7. The measurement agrees with the LHCb result. All 
ratios are well described by the recent perturbative QCD predictions [8]. The expectations from models 
in Refs. [3, 5, 7] as well as the sum-rules prediction [4] for the ratio 'Rd*+id+ consistent with the 
measurement. The QCD relativistic potential model predictions [3] are consistent with the measured 
'Rd+/ 7 t* while the expectations from the sum rules [4] and models in Refs. [5-7] are somewhat 
smaller than the measured value. The predictions in Refs. [3-5, 7] are also generally smaller than the 
measured ratio however, the discrepancies do not exceed two standard deviations when taking 

into account only the experimental uncertainty. 


Table 5; Comparison of the results of this measurement with those of LHCb [10] and theoretical predictions based 
on a QCD relativistic potential model [3], QCD sum rules [4], relativistic constituent quark model (RCQM) [5], 
BSW relativistic quark model (with fixed average transverse quark momentum to - 0.40 GeV) [6], light-front 
quark model (LFQM) [7], perturbative QCD (pQCD) [8], and relativistic independent quark model (RIQM) [9]. 
The uncertainties of the theoretical predictions are shown if they are explicitly quoted in the corresponding papers. 
Statistical and systematic uncertainties added in quadrature are shown for the results of ATLAS and LHCb. 





r±±/F 

Ref. 

3.8 ± 1.2 

2.90 ± 0.62 

10.4 ± 3.5 

1 8+1-2 
2-8-09 

2.37 ± 0.57 

0.38 ± 0.24 

0.52 ± 0.20 

ATLAS 

LHCb [10] 

2.6 

1.3 

2.0 

2.2 

2.06 ± 0.86 

O ^^+0.49 

4.5 

5.2 

5.7 

1.7 

3.9 

2.9 

3.01 ± 1.23 

p 54+0.07 
-0.21 

0.48 ± 0.04 

0.410 

QCD potential model [3] 
QCD sum rules [4] 
RCQM [5] 
BSW [6] 
LFQM [7] 
pQCD [8] 
RIQM [9] 


The measured fraction of the A++ component agrees well with the prediction of the relativistic independ¬ 
ent quark model [9] and perturbative QCD [8] . 

9 Conclusion 

A study of and decays has been performed. The ratios of the branching 

fractions ^he transverse po¬ 

larisation fraction of B^ —> JI decay have been measured by the ATLAS experiment at the LHC 
using pp collision data corresponding to an integrated luminosity of 4.9 fb“^ at 7 TeV centre-of-mass 
energy and 20.6 fb“^ at 8 TeV. The polarisation is found to be well described by the available theoretical 
approaches. The measured ratios of the branching fraction are generally described by perturbative QCD, 
sum rules, and relativistic quark models. There is an indication of underestimation of the decay rates 
for the B^ decays by some models, although the discrepancies do not exceed two standard 

deviations when taking into account only the experimental uncertainty. The measurement results agree 
with those published by the LHCb experiment. 
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Figure 7: Comparison of the results of this measurement with those of LHCb [10] and theoretical predictions based 
on a QCD relativistic potential model [3], QCD sum rules [4], relativistic constituent quark model (RCQM) [5], 
BSW relativistic quark model (with fixed average transverse quark momentum o) - 0.40 GeV) [6], light-front 
quark model (LFQM) [7], perturbative QCD (pQCD) [8], and relativistic independent quark model (RIQM) [9]. 
The uncertainties of the theoretical predictions are shown if they are explicitly quoted in the corresponding papers. 
Statistical and systematic uncertainties added in quadrature are quoted for the results of ATLAS and LHCb. 
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